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On the basis of the perturbation theory in the fine structure constant a and the 
ratio of the electron to muon masses we calculate one-loop vacuum polarization and 
electron vertex corrections and the nuclear structure corrections to the hyperfine 
splitting of the ground state of muonic helium atom (^e^-ffe). We obtain total result 
for the ground state hyperfine splitting l^v^f'^ = 4465.526 MHz which improves the 
previous calculation of Lakdawala and Mohr due to the account of new corrections. 
The remaining difference between the theoretical result and experimental value of 
the hyperfine splitting equal to 0.522 MHz lies in the range of theoretical error and 
requires the subsequent investigation of higher order corrections. 

PACS numbers: 31.30.Jv, 12.20.Ds. 32.10.Fn 



I. INTRODUCTION 

Muonic helium atom {^e^He) represents the simple three-body atomic system. The 
interaction between magnetic moments of the muon and electron leads to the hyperfine 
structure (HPS) of the energy levels. The investigation of the energy spectrum of this 
three-particle bound state is important for the further check of quantum electrodynamics. 
Hyperfine splitting of the ground state of muonic helium was measured many years ago with 



sufficiently high accuracy 



groun( 
18: 



exp 



4465.004(29) MHz. 



(1) 



Contrary to the energy levels of two-particle bound states which were accurately calculated 
in quantum electrodynamics [l, 0, [s], 0, 0, B 0], the hyperfine splitting of the ground state 
in muonic helium atom was calculated on the basis of the perturbation theory (PT) and the 
variational method with essentially less accuracy 



Indeed, the theoretical errors of the results obtained in Refs. [10, M, M, M, M, M, M, M, M, 
19|, |20|, |21j lie in the interval 0.05-7-1.8 MHz. The variational method gives the high numerical 



accuracy of the calculation what was demonstrated in Refs. [15|, Il8|, l20|, |21[. But higher order 



corrections were accounted in this approach less precisely. So, for instance, the theoretical 
uncertainty 0.05 MHz in Ref. 1^ is estimated only from the numerical convergence of the 
results obtained by the variational method for the lowest order hyperfine splitting. The 
nonvariational calculation of the lowest order contribution to the hyperfine structure was 
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performed in the hyperspherical harmonic method in Ref. 22(]. But numerous important 
corrections to the hyperfine sphtting which are necessary for the successful comparison with 
the experimental data were not considered in 2^ . 

Many theoretical efforts were focused on the calculation of different corrections to the 
Fermi energy which is of the fourth order over the fine structure constant a. The first part of 
the calculations was devoted to the recoil corrections which contain the ratio of the electron 
and muon masses [13, 17|. The second group was connected with the relativistic and QED 
effects which include another small parameter a IOI. IIII. 14 The particle coordinates in the 



muonic helium atom were taken in Refs. 0, OTBrilUS, H H E, Q by the various 
ways. So, numerical results for the corrections of different order over the fine structure 
constant a and the ratio of the particle masses obtained in these papers are difficult for the 
direct comparison. 

The bound particles in muonic helium atom have different masses me <^ <^ mo,. 
As a result the muon and a-particle compose the pseudonucleus (/ig-ffe)"^ and the muonic 
helium atom looks as a two-particle system in the first approximation. The perturbation 
theory approach to the investigation of the hyperfine structure of muonic helium based on 
the nonrelativistic Schrodinger equation was developed previously by Lakdawala and Mohr 
in Refs. fl^. Three-particle bound system (/iCg-ffe) is described by the Hamiltonian: 



Hn 
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where and Xg are the coordinates of the muon and electron relative to the helium nucleus. 
Me = rrif^ma/ {jrif. + ^ah.Mn = '^n^jna/ {m^ + rUa) are the reduced masses of subsystems 
{e\He)^ and {ii\He)^ 12|, ll7|]. The hyperfine part of the Hamiltonian is 



5(x^ - Xe 



(4) 



where cr^ and cr^ are the spin matrices of the electron and muon, Ke and are the electron 
and muon anomalous magnetic moments. In the initial approximation the wave function of 
the ground state has the form [H, 17]: 



71 



■aMeXe 



(5) 



Then the basic contribution to the singlet-triplet hyperfine splitting can be calculated ana- 
lytically from the contact interaction (4): 
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Numerically the Fermi splitting is equal vp = 4516.915 MHz. We express further the hy- 
perfine splitting contributions in the frequency unit using the relation AE^^" = 2TihAv^^^ . 
The recoil correction determined by the ratio Me/M^ in Eq.(6) amounts Av^l^ = —33.525 
MHz [l2|. Modern numerical values of fundamental physical constants are taken from 
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the 



the electron mass mg = 0.510998918(44) 



paper 

0.1056583692(94) GeV, the fine structure constant a~ 



10 ^ GeV, the muon mass 
137.03599911(46), the he- 



hum mass m{\He) = 3.72737904(15) GeV, the electron anomalous magnetic moment 
1.1596521869(41) ■ 10"'^, the muon anomalous magnetic moment = 1.16591981(62) ■ lO"'^. 

Analytical and numerical calculation of corrections which are determined by the Hamilto- 
nians AH and Aif^ec in the second order perturbation theory was performed in Refs. jl2, 17 . 
Their results and the order of the calculated contributions are presented in Table I. In this 
work we aim to refine the calculation of Lakdawala and Mohr using their approach to the de- 
scription of the muonic helium atom. A feature that distinguishes light muonic atoms among 
the simplest atoms is that the structure of their energy levels depends strongly on the vac- 
uum polarization, nuclear structure and recoil effects jl, 0, H, B, 0, B [9j. So, we investigate 
such contributions of the one-loop electron vacuum polarization of order a^Me/M^ and the 
nuclear structure of order which are significant for the improvement of the theoretical 
value of the hyperfine splitting. Another purpose of our study consists in the improved cal- 
culation of the electron one-loop vertex corrections to HFS of order using the analytical 
expressions of the Dirac and Pauli form factors of the electron. 





FIG. 1: The vacuum polarization effects. The dashed line represents the Coulomb photon. The 
wave line represents the hyperfine part of the Breit potential. G is the reduced Coulomb Green's 
function. 



II. EFFECTS OF THE VACUUM POLARIZATION 

The vacuum effects change the interaction (2)- (3) between particles in muonic helium 
atom. One of the most important contributions is determined by the one-loop vacuum 
polarization (VP) and electron vertex operator. The electron vacuum polarization and 
vertex corrections to the hyperfine splitting of the ground state contain the parameter equal 
to the ratio of the Compton wave length of the electron and the radius of the Bohr orbit in the 
subsystem {^^He)^'- ni^a/mf.= 1.50886. . . . The expansion over a for such contributions 
to the energy spectrum is unfit to use. So, we calculate them performing the analytical 
or numerical integration over the particle coordinates and other parameters without an 
expansion in a. The effect of the electron vacuum polarization leads to the appearance of a 



4 



number of additional terms in the interaction operator which we present in the form [s], [13] : 
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where = |xe — x^|. They give contributions to the hyperfine sphtting in the second 
order perturbation theory and are discussed below. In the first order perturbation theory 
the contribution of the vacuum polarization is connected with the modification of the hy- 
perfine splitting part of the Hamiltonian (4) (the diagram (a) in Fig. 1). In the coordinate 
representation it is determined by the integral expression 25, [1^, 27]: 



STJiem^ 4 Stt 
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Averaging the potential (10) over the wave function (5) we obtain the following contribution 
to the hyperfine splitting: 
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There are two integrals over the muon and electron coordinates in Eq.(ll) which can be 
calculated analytically: 



h = j dxe J dXf,e~ 
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They are divergent separately in the subsequent integration over the parameter ^. But their 
sum is finite and can be written in the integral form: 



p{m- 



Me _|_ o rrteg Me _j_ meg / 9 _|_ meg 
2M^ 2AI^a2M^ ~'~ 2M^a \ ' 2M^o 



f 1 _|_ meg A j\/e _|_ meg 

V"^ 2Mf.ay \2M^ 2M^a 



0.035 Mi/z. 



(14) 



The order of this contribution is determined by two small parameters a and M(,/M^ which 
are written explicitly. The correction Az/^p is of the fifth order in a and the first order in the 
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ratio of the electron and muon masses. The contribution of the muon vacuum polarization to 
the hyperfine splitting is extremely small (~ 10~^ MHz). One should expect that two- loop 
vacuum polarization contributions to the hyperfine structure are suppressed relative to the 
one-loop VP contribution by the factor a/vr. This means that at present level of accuracy 
we can neglect these corrections because their numerical value is not exceeding 0.001 MHz. 
Higher orders of the perturbation theory which contain one-loop vacuum polarization and 
the Coulomb interaction (3) lead to the recoil corrections of order PpOi^ In Such terms 
which can contribute 0.004 MHz are included in the theoretical error. 

It is useful to compare the obtained result (14) with the calculation of VP contribution 
to the HFS in which the expansion in a is used. Instead of the potential (10) we obtain the 
following operator in the coordinate representation: 

AVV'^''(x,J = — ^!^^i-v25(xeJ. (15) 

Then the contribution of (15) to the hyperfine structure can be derived in the analytical 
form: 

~hfs_ 8a (aM^\^ M, 1 



Az>^^/ = uf—[ ^ — 1 = 0.060 MHz. (16) 

157r V me / M^fij^ m^A^ ^ ^ 

This calculation demonstrates the need to employ the exact potentials (7)-(9) for the study 
of the electron vacuum polarization corrections. 

Let us consider corrections of the electron vacuum polarization (7)-(9) in the second 
order perturbation theory (SOPT) (the diagram (b) in Fig.l). The contribution of the 
electron-nucleus interaction (7) to the hyperfine splitting can be written as follows: 

soft ea = 3^^^ / ^^^i j rfxa j ^^^^ pi.^) di%o{^^Weo{^^) X (17) 



X 



^^„(x3)V^en'(x3)Cn(x2)^e*n'(xi) 



Here the summation is carried out over the complete system of the eigenstates of the electron 
and muon excluding the state with n, n' = 0. The computation of the expression (17) is 
simplified with the use of the orthogonality condition for the muon wave functions: 



X 



(18) 






Me 



- In ( ]+Ei[ ] +--2C- t^{xi + X3) + 



2M/ 



0.151 MHz, 



where a;< = min(xi,X3), a;> = max(a;i, xs), C = 0.577216... is the Euler's constant and 
Ei{x) is the exponential- integral function. It is necessary to emphasize that the transforma- 
tion of the expression (17) into (18) is carried out with the use of the compact representation 



for the electron reduced Coulomb Green's function obtained in Refs. |12| . 128 

^en(x3)V^:„(Xi) _ ^^^(^^+^3) 



Ge(xi, X3) = J2 p ' p' ' = — -e 



2aMeX^ 



(19) 



6 



7 1 — g2aMe2;< 

- ln(2aMeX>) - ln(2aMeX<) + Eii2aM^x^) + - - 2C - aMeixi + 0:3) + 

2 2aMpX 



The contribution (18) has the same order of the magnitude 0{a^^) as the previous cor- 
rection (14) in the first order perturbation theory. Similar calculation can be performed in 
the case of muon-nucleus vacuum polarization operator (8). The intermediate electron state 
is the IS state and the reduced Coulomb Green's function of the system appearing in the 
second order PT transforms to the Green's function of the muon. The correction of the 
operator (8) to the hyperfine splitting is obtained in the following integral form: 



^^vp SOFT ^a = ^f:^ piOdU x^dxs X2dx2e V We V 
^ Stt Ji Jo Jo 



(20) 
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0.048 MHz. 



The vacuum polarization correction to HPS which is determined by the operator (9) in the 
second order perturbation theory is the most difficult for the calculation. Indeed, in this 
case we have to consider the intermediate excited states both for the muon and electron. 
PoUowing Ref. [13] we have divided total contribution into two parts. The first part in which 
the intermediate muon is in the IS state can be written as: 

xj^ xldxie-''^^'^+^^'-^^' p{Od^AVvF^{xi)Ge{xi,xs), 
where the function Vyp fi{xi) is equal 

AVvF u{Xl) = f rfx^e-^"*^^-^ (2«^m)' « g-2m.g|x,-x,| ^ (22) 

J n |xi — X2I 

After the substitution (22) in (21) the numerical integration gives the result: 

soft f^ein = 0) = -0.030 MHz. (23) 

Second part of the vacuum polarization correction to the hyperfine splitting due to the 
electron-muon interaction (9) can be presented as follows: 



A4^; soft ,e{n ^0) = J d^-s J d^2 f p{Od^r,oi^,)reoi^,) X (24) 



16a 

A/f p-^'|x3-xi| ^ 

X E ^.n(x3)^;„(x2)-^ I II 1 e-^-^«l-^-l^,0(x2)^e0(Xl). 

n^O 1X3 — Xi I |X2 — Xi| 

In the expression (24) we have replaced the exact electron Coulomb Green's function by 
the free electron Green's function which contains b = [2Me(i?^„ — -E^o ~ E^oY^'^. (see more 
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detailed discussion of this approximation in Re fs.[l2l . [l7| ) . We also replace the electron wave 
functions by their values at the origin as in Ref. [12] neglecting higher order recoil corrections. 
After that the integration over xi can be done analytically: 
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(ixi 



g-fe|x3-Xl| g-2TOeS|x2-Xl| 

1X3 -Xi I 1X2 -Xi I 



47r 



1 
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\ _ g-2mcC|x3-X2| 



|x3 — X2I 6^ — Aml^"^ 

-2me5|x3-X2h ^2 



-2me5|x3-X2 



+ 



(25) 



63 63(^3 _xi)2 



&^|X3 - X2I 



+ ... 



where we have performed the expansion of the first exponential in the square brackets over 
powers of 6|x3 — X2I. As discussed in Ref. 12| one can treat this series as an expansion over 



the recoil parameter yMe/M^. For the further transformation the completeness condition 
is useful: 

53 ?/'^„(x3)V';^„(x2) = (5(X3 - X2) - ^^o(x3)^*o(x2)- (26) 

The wave function orthogonality leads to the zero results for the second and fifth terms in 
the square brackets of Eq.(25). The first term in Eq.(25) gives the leading order contribution 
in two small parameters a and Me/M^: 



^^VP SOFT ue(^ 7^ 0) = Az/ii + Az/12, Az/n 



-Up, 



(27) 



Al/i2 = Up 



2a' 
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e (4«M^ + 2m,0' 



256 



rripC mli'^ mli^ 

232^^^ + 80—1^ + 10 



(28) 

The numerical value of the sum Auu+Auu is included in Table I. It is important to calculate 
also the contributions of other terms of the expression (25) to the hyperfine splitting. Taking 
the fourth term in Eq.(25) which is proportional to 6^ = 2Me(-E^„ — E^q) we have performed 
the sequence of the transformations in the coordinate representation: 



/ dX2 / C/X3?/;*0(X2)V^^„(X3)?/;*„(X2)|X3 - X2|^^0(X2) 

n=0 



(29) 



J dX2 J (iX35(x3 - X2) 
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|X3 -X2|^*0(X3 



^mo(x2). 



Evidently, we have the divergent expression in Eq.(29) due to the presence of the ^-function. 
The same divergence occurs in the other term containing b' in the square brackets of Eq.(25). 
But their sum is finite and can be calculated analytically with the result: 



Up 



a'' Ml' 
5 



(30) 
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Numerical value of this correction 0.0002 MHz is essentially smaller than the leading order 
term. Let us consider also the nonzero term in Eq.(25) proportional to b^. First of all, it 
can be transformed to the following expression after the integration over ^: 



-Up 



457r' 



fMlMl 
M„ ml 



where the sum 5*3/2 is defined as follows: 



(31) 



(32) 



l/2aM^. Using the known analytical expressions for the dipole matrix 



elements entering in Eq.(32) in the case of the discrete and continuous spectrum |3l. |29 
can write separately their contributions to the sum 5*3/2 in the form: 

°^ 28n^(n - 1)2"-- 
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(33) 



(34) 



As a result 53/2 = 3.2722 .... The similar calculation of the sum 5i/2 relating to this 
problem (see Ref. [13]) gives 5i/2 = 2.9380 . . .. Numerical value (31) is taken into account in 
the total result presented in Table I. 




FIG. 2: Vacuum polarization effects in the second order perturbation theory. The dashed line 
represents the first part of the potential AH (3). The wave line represents the hyperfine part of 
the Breit potential. 



There exists another contribution of the second order perturbation theory in which we 
have the vacuum polarization perturbation connected with the hyperfine splitting part of the 
Breit potential (10) (see Fig.2). Another perturbation potential in this case is determined by 
the first term of relation (3). We can divide this correction into two parts as previously. One 
part with n = corresponds to the ground state muon. The other part with n ^ accounts 
the excited muon states. The 5-function term in Eq.(lO) gives the following contribution to 
HFS at n = (compare with Ref. [l^): 



^^V^P SOFT lli^ 



0) = / 



(35) 
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Obviously, this integral in the variable ^ is divergent. So, we have to consider the contribution 
of the second term of the potential (10) to the hyperfine splitting which is determined by 
the more complicated expression: 

a4^; soft i2{n = 0) = ^ / piOed^ / rfx3^eo(x3)AV-i(x3) X (36) 

^ ^en'(x3)^eV(xi) ^ / \l ( \ 
^ 2^ AV2(Xl)^e0(xl), 



where 



Ayi(x3) = / dX4V^;o(x4) ^ ° ' I ^M0(X4) = (37) 
J 1X3 — X4I 



X3[(4«M^)2-(2meO']' 
Ay2(xi) = / rfx2^^o(x2) f i " I - -\ ^,o(x2) = "-(1 + 2aM ^x,)e''-^ (38) 

J \ Xo — Xi Xi / Xi 



Nevertheless, integrating over all coordinates in Eq(36) we obtain the following result in the 
leading order with respect to the ratio (Mg/M^): 

A hfs , r,. m,Ml ,^^^^^ 32 + 637 + 447^ + 117^ 
AVp soft An = 0) = i^F^^r^-rpi L PiO^^^ T^TTXI ' (39) 



Me 967rM2 Ji "^'^'^ ^ (1 + 7)^ 

where 7 = meC,/'2aM^. This integral also has the divergence at large values of the parameter 
^. But the sum of integrals (35) and (39) is finite: 

soft 11 = 0) + Auy-^p SOFT 12 = 0) = 0.008 MHz. (40) 



G 



FIG. 3: Nuclear structure effects in the second order perturbation theory. The bold point represents 
the nuclear vertex operator. The wave line represents the hyperfine part of the Breit potential. 



Let us consider now the terms with n ^ 0. The delta-like term of the potential (10) gives 



the contribution to HPS known from the calculation of Ref. 12 



^^VF SOFT 21 ('^ 7^ 0) 
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35Me 
16M„, 



(41) 
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Another correction from the second term of the expression (10) can be simphfied after the 
replacement the exact electron Green's function by the free electron Green's function: 



A4^; soft 22{n 7^ 0) = - 'gJJ^^- f piOed^ J J d^-, x (42) 

/g-2me,^|x3-x4| oo 
d^4i^*^,oi^4)—, 7— V'A<n(x4)V'Mn(x2)|x3 - X2|V^^o(x2) 
Xa X4 „ /„ 



The analytical integration in Eq.(42) over all coordinates leads to the result: 

aM, 



^^vp SOFT 22(^ ^ 0) = -^F-^;^ ^ p{Odi 



1 1 1 _ 21572 357^ 

'4+- + IO7 + 



7 (1 + 7)4 \^ 7 ' 16 16 ^ 

(43) 

The sum of expressions (41) and (43) gives again the finite contribution to the hyperfine 
splitting: 

^^VP SOFT 2lin 7^ 0) + Az/Jp SOFT 22(^ 7^ 0) = (44) 

37iM^Ji ^ 16(1 + 7)4 

Despite the fact that the absolute values of the calculated VP corrections (23), (27), (28), 
(30), (31), (40), (44) are sufficiently large , their summary contribution to the hyperfine 
splitting (see Table I) is small because they have different signs. 



III. NUCLEAR STRUCTURE AND RECOIL EFFECTS 

Another significant corrections to the hyperfine splitting of muonic helium atom which 
we study in this work are determined by the nuclear structure effects. They are specific for 
any muonic atom. In the leading order over a they are described by the charge radius of 
a-particle r^. If we consider the interaction between the muon and the nucleus then the 
nuclear structure correction to the interaction operator has the form (^: 

AVstrA^,) = '^TcZa <rl> 5{v^). (45) 

The contribution of the operator ^Vstr,^l to the hyperfine splitting appears in the second 
order perturbation theory (see the diagram in Fig. 3). In the beginning we can write it in 
the integral form: 

^^1:, = ^^rl^ {2aM,f' j d^,r,oi?^,)\^eo{^,)\'G,{-K,, 0, E,,). (46) 

After that the analytical integration over the coordinate X3 in Eq.(46) can be carried out 
using the representation of the muon Green's function similar to expression (17). The result 
of the integration of order 0{a^) is written as an expansion in the ratio Mf,/M^: 

Az^^t, = -i^rla^Mlrl " + • • •) = -0-007 MHz. (47) 
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TABLE I: Hyperfine singlet-triplet splitting of the ground state in the muonic helium atom. 



Contribution to the HFS 


Au^J^', MHz 


Reference 


The Fermi splitting 


4516.915 


(6), [12j 


Recoil correction of order 


-33.525 


(6),112J 


Correction of muon anomalous 
magnetic moment of order ot' 


/I /I 

0.z44 




Recoil correction of order 

l\/r Im \ 1 ( A/T I'h/T \ 

a {ivie/ma)y [Mg/M^) 


0.079 


[12, la 


Correction due to the perturbation 
(3) in the second order PT 
of order o;^^^ 


-29.650 


[12,17] 


Relativistic correction of order 


0.040 


[10] 


One-loop VP contribution in I7 
(e/i)-interaction of order a^^^ 


0.035 


(14) 


One-loop VP contribution in the 
electron-muon interaction in the 
second order PT of order a'^jj^ 


-0.145 


(23),(27),(28),(30), 
(31),(40),(44) 


One-loop VP contribution in the 
electron-nucleus interaction in the 
second order PT of order a^^^ 


0.151 


(18) 


One-loop VP contribution in the 
muon-nucleus interaction in the 
second order PT of order a^jj^ 


0.048 


(20) 


Nuclear structure correction of order 
Ct in Liie oeconcL oiciei r± 


-0.010 


(47),(49) 


Recoil correction of order 
a^ime/nifj,) ln(me/m^) 


0.812 


(51) m 


Vertex correction of order 


-0.606 


[5, 31,32,33] 


Electron vertex contribution 
of order 


6.138 


(53),(56),(58), 
(59),(61),(62) 


Summary contribution 


4465.526 





Numerical value of the contribution AuJ/^^ is obtained by means of the charge radius of 
the a-particle = 1.676 fm. The same approach can be used in the study of the electron- 
nucleus interaction. The electron feels as well the distribution of the electric charge of a 
particle. The corresponding contribution of the nuclear structure effect to the hyperfine 
splitting is determined by the expression: 

^^str,e = rl j dXi j rfX3|^/'*o(x3)PV^eo(x3)G'^(x3,Xi,£'eo)^/'eo(Xl)5(Xi). (48) 
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Performing the analytical integration in Eq.(46) we obtain the following series: 



Az/: 



hfs 
str,e 



4 2 71^-2 2 



5 - In ^ 

M, Ml 



3 In 



Me 



-7 + 



Ml 17 
Ml [2 



- 31n 



Me 

M,, 



(49) 



-0.003 Mi/z. 



We have included in Table I the total nuclear structure contribution which is equal to the 
sum of the numerical values (45) and (47). 



a b 
FIG. 4: Two photon exchange amplitudes in the electron-muon hyperfine interaction. 

Special attention has to be given to the recoil corrections connected with two-photon 
exchange diagrams shown in Fig.4 in the case of the electron-muon interaction. For the 
singlet-triplet splitting the leading order recoil contribution to the interaction operator be- 
tween the muon and electron is determined as follows 



9 

AKei%(x,e) = 8^ In ^5(x,e). (50) 

— me 

Averaging the potential AVjll^^e o'^^r the wave functions (5) we obtain the leading order 
recoil correction to the hyperfine splitting: 

. hf^ 3q; mpm., , m,, 

Az/^Yc = In — = 0.812 MHz. (51 

IT — me 

There exist also the two-photon interactions between the bound particles of muonic helium 
atom when one hyperfine photon transfers the interaction from the electron to muon and 
another Coulomb photon from the electron to the nucleus (or from the muon to the nucleus). 
Supposing that these amplitudes give smaller contribution to the hyperfine splitting we 
included them in the theoretical error. 



IV. ELECTRON VERTEX CORRECTIONS 



In the initial approximation the potential of the hyperfine splitting is determined by 
Eq.(4). It leads to the energy splitting of order a'^. In QED perturbation theory there is 
the electron vertex correction to the potential (4) which is defined by the diagram in Fig. 5 
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(a). In momentum representation the corresponding operator of hyperfine interaction has 
the form: 



6a 



(T (T 

e /i 



(52) 



where G^^{k'^) is the electron magnetic form factor. We extracted for the convenience the 



factor a/vr from G^^{k 



factor G^^){k'^ 



M 

Gl^i(O) 



Usually used approximation for the electron magnetic form 



My^ ) vj; — 1 + is not quite correct in this task. Indeed, characteristic 

momentum of the exchanged photon is /c ~ aM^. It is impossible to neglect it in the 
magnetic form factor as compared with the electron mass nie- So, we should use exact 



one-loop expression for the magnetic form factor which was obtained by many authors [24 
Let us note that the Dirac form factor of the electron is dependent on the parameter of the 
infrared cutoff A. We take it in the form A = rriea using the prescription rrieOp' -C A -C rrie 
from Ref. 



G 



G 



a 



FIG. 5: The electron vertex corrections. The dashed Une represents the Coulomb photon. The 
wave line represents the hyperfine part of the Breit potential. G is the reduced Coulomb Green's 
function. 



Using the Fourier transform of the potential (52) and averaging the obtained expression 
over wave functions (5) we represent the electron vertex correction to the hyperfine splitting 
as follows: 



a 



nif. 



vertex 



Up 



327r2 

2 



k^dk 



k' 



' Me 

2M„ 



+ 



/ me 



1 



(53) 



4.214 MHz. 



Let us remark that the contribution (53) is of order a^. Numerical value (53) is obtained 
after numerical integration with the one-loop expression of the electron magnetic form factor 
^m(^^)- If ^6 use the value G^^}{k'^ = 0) then the electron vertex correction is equal 
5.244 MHz. So, using the exact expression of the electron form factors in the one- loop 
approximation we observe the 1 MHz decrease of the vertex correction to the hyperfine 
splitting from I7 interaction. Taking the expression (52) as an additional perturbation 
potential we have to calculate its contribution to HFS in the second order perturbation 
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theory (see the diagram in Fig. 5(b)). In this case the dashed hne represents the Coulomb 
Hamiltonian AH (3). Following the method of the calculation formulated in previous section 
(see also Refs.[ll, [l3|) we divide again total contribution from the amplitude in Fig.5(b) into 
two parts which correspond to the muon ground state {n = 0) and muon excited intermediate 
states {n ^ 0). In this way the first contribution with n = takes the form: 



^J^virtexin = 0) 



6a 



dk I c/xi / (iX3?/^eo(x3) X (54) 



xAVi{k, X3)Ge(xi, X3)AV2(xi)^eo(xi) 

where AV2(xi) is defined by Eq.(38) and 



Ari(A;,X3) 



sin(A;|x3 -X4I) 

dX4?/'^0 (X4) , , ^^0 (X4 J 



sm 



X3 - X4 



X2. 



1 + 



fc2 



(4aM^)2 



(55) 



Substituting the electron Green's function (19) in Eq.(54) we transform desired relation to 
the integral form: 



^^vlrtexin = 0) = Up 



a 



167r2 \aM„ 



Me 
M„ 



k 



Gg(fc^)-1 



dk 



1 + 



(4aM^)2 



X 



(56) 



X 



roo _JW^ . / rriek \ Xx\ ~xx{\+^\ , 



X 



2M^ 



- ln(-— — x<) - ln(-— — a;>) + Ei[-^^^x<) + - - 26 - tt^I^^i + 2^3) + 



■a;< 



'2M„ 



^2M^ 
-0.210 Mi/z. 



4M„ 



One integration over the coordinate Xi is carried out analytically and two other integrations 
are performed numerically. Second part of the vertex contribution (Fig. 5(b)) with n 7^ can 
be reduced to the following form after several simplifications which are discussed in section 
II (see also Refs.fii, Q): 



hfs Sa^M^Mf^ 

^^vLtexin ^ 0) = Up — ^ 



-2aMnX2 



X 



|^"fcsin(A;|x3-X4|)(Gi^^(fc2)-l 



(iX2 



|X3 - X2I 



-aMeX3 



rfxs 



-2aMuX4 



C?X4 X 



(57) 



[5(X4 - X2) - V^^o(x4)V^^.o(x2)] • 



1X3 - X4I 

We divide expression (57) into two parts as provided by two terms in the square brackets of 
(57). After that the integration (57) over the coordinates xi, X3 is carried out analytically. 
In the issue we obtain (71 = Me/4M^, 72 = rriek/AaM^): 



^ilrtexi.'^ 7^ 0) — 



327r2 \aMf,J Jo 



1 



dk X 



(58) 



471(7? - 1) 7i(3 + 7?) , 47?(7? - 1) ^ 1 + 37? 



2\3 



(I + 7IP (7? + 72T (7? + 7l)^ 



2.516 MHz, 
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^^2!vertex{n 7^ 0) = -Up 



a 



rrir. 



2\2 



327r2 
(71 + 1) 



dk X 



(59) 



(7? + 7l) [(l + 7i)^ + 7iF (71 + 1)^ + 71 ilf + lir 



-0.831 Mi/^. 



It is necessary to emphasize that the theoretical error in the contributions Aui 



hfs 



2,vertex 



[n ^ 0) 

is determined by the factor ^ M(,/M^ connected with the omitted terms of the expansion 
similar to Eq.(25) (see also Refs. [l^. 17|). It can amount to 10% of the results (58), (59) 
that is the value near 0.2 MHz. 

Until now we consider the electron vertex corrections connected with the hyperfine part 
of the interaction Hamiltonian (4). But in the second order perturbation theory we should 
analyze vertex corrections to the Coulomb interactions of the electron and muon, electron 
and nucleus. Then in the coordinate representation we have the following potential: 



A T/*^ ( 

vertex, eN > 



~^ ^^vertex,efj.{-^eij) 



2a' 



Gg^(A:2 



k 



■dk 



( sin(fca;, 



- 2 



sin(fcxe 



X. 



(60) 

where we extract again the factor a/ir from G^^\k'^) — 1 . G^'* is the electron electric 
form factor. One part of the contribution in Fig. 5(c) is specified by the electron-muon 
intermediate states in which the muon is in the ground state n = 0. This correction is 
determined by both terms in the round brackets of Eq.(60) and can be presented as follows: 




= -1.321 MHz. 

The index "C" means that the vertex correction to the Coulomb part of the Hamiltonian 
is considered. Excited states of the muon (n ^ 0) contribute to another part of the matrix 
element (Fig.5(c)). Changing the Coulomb Green's function of the electron by free Green's 
function (see discussion in section II) we can make the coordinate integration and express 
the correction to HFS as one-dimensional integral: 



^^C.vertexi^ 7^ 0) = -Up 



8a /aM, 



G^^\k^ 



dk 



1 + 



(62) 



1.770 MHz. 



The electron vertex corrections investigated in this section have the order a^ in the hyperfine 
interval. Summary value of all obtained contributions (53), (56), (58), (59), (61), (62) is 
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equal to 6.138 MHz. It differs by a significant value 0.894 MHz from the result 5.244 MHz 
which was used previously by many authors for the estimation of the electron anomalous 
magnetic moment contribution. 



V. CONCLUSIONS 

In the present study, we have performed the analytical and numerical calculation of 
several important contributions to the hyperfine splitting of the ground state in muonic 
helium atom connected with the vacuum polarization, the nuclear structure effects and 
the electron vertex corrections. To solve this task we use the method of the perturbation 
theory which was formulated previously for the description of the muonic helium hyperfine 
splitting in Refs. l3, 17 1. We have considered corrections of order oc^ of the electron vacuum 



polarization and electromagnetic form factors and the nuclear structure effects of order 
a^. The numerical values of the corresponding contributions are displayed in Table I. We 
present in Table I the references to the calculations of other corrections which are not 
considered here. The relativistic correction was obtained in Ref. 10 , the vertex correction 



was calculated in the case of hydrogenic atoms in Refs.js], [3l|, 32, 33|. Basic contributions 
to the hyperfine splitting obtained by Lakdawala and Mohr are also included in Table I 
because our calculation is closely related to their approach. 
Let us list a number of features of the calculation. 

1. For muonic helium atom, the vacuum polarization effects are important and give 
rise to the modification of the two-particle interaction potential which provides the oc^^- 
order corrections to the hyperfine structure. The next to leading order vacuum polarization 
corrections (two-loop vacuum polarization) are negligible. 

2. The electron vertex corrections should be considered with the exact account of the 
one-loop electromagnetic form factors of the electron because the characteristic momentum 
incoming in the electron vertex operator is of order of the electron mass. 

3. In the a^-order the nuclear structure corrections to the ground state hyperfine splitting 
are expressed in terms of the charge radius of a-particle. 

4. Analyzing the one-loop electron vacuum polarization and vertex effects and the nu- 
clear structure contributions in each order in a, we have taken into account recoil terms 
proportional to the ratio of the electron and muon masses. 

The resulting numerical value 4465.526 MHz of the ground state hyperfine splitting in 
muonic helium is presented in Table I. It is sufficiently close both to the experimental result 
(1) and the earlier performed calculations by the perturbation theory, variational approach 
and the Born-Oppenheimer theory: 4464.3±1.8 MHz [l75, 4465.0±0.3 MHz [Te^, 4462.9 MHz 
Q, 4450.4 ±0.4 MHz Q, 4459.9 MHz ^], 4464.87 ±0.05 MHz |19]. The estimation of the 
theoretical uncertainty can be done in terms of the Fermi energy vp and small parameters 
a and the ratio of the particle masses. On our opinion there exist several main sources 
of the theoretical errors. First of all, as we mentioned above comprehensive analytical 

and numerical calculation of recoil corrections of orders o^tt^, a^T?!, ^^"^ ln(M„/Me) was 

— I 

carried out by Lakdawala and Mohr in the second order PT in Refs. |l2l . Il7l | . The error of 
their calculation connected with the correction z/^|pln^ consists 0.6 MHz. The second 

source of the error is related to contributions of order o?vp 0.2 MHz which appear 
both from QED amplitudes and in higher orders of the perturbation theory. Another part 
of the theoretical error is determined by the two-photon three-body exchange amplitudes 
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mentioned above. They are of the fifth order over a and contain the recoil parameter 
{nif^/ma) ln(me/ma), so that their possible numerical value can be equal ±0.05 MHz. Finally, 
a part of theoretical error is connected with our calculation of the electron vertex corrections 
of order in section IV. It consists at least 0.2 MHz (see the discussion after Eq.(59)) We 
neglect also the electron vertex contributions of order Vi,-aMe/M^ k, 0.2 MHz which appear 
in higher orders of the perturbation theory. Thereby, the total theoretical uncertainty is 
not exceeded ±0.7 MHz. The existing difference between the obtained theoretical result 
and experimental value of the hyperfine splitting (1) equal to 0.522 MHz hes in the range 
of total error. Theoretical error which remains sufficiently large in the comparison with 
the experimental uncertainty, initiates further theoretical investigation of the higher order 
contributions including more careful construction of the three-particle interaction operator 
connected with the multiphoton exchanges. 
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